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Interacting Electron Wave Packet Dynamics in a Two-dimensional Nanochannel 
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Classical and quantum dynamics are important limits for the understanding of the transport 
characteristics of interacting electrons in nanodevices. Here we apply an intermediate semiclassical 
approach to investigate the dynamics of two interacting electrons in a planar nanochannel as a 
function of Coulomb repulsion and electric field. We find that charge is mostly redistributed to 
the channel edges and that an electric field enhances the particle-like character of electrons. These 
results may have significant implications for the design and study of future nanodevices. 



Introduction — Ideal ultrasmall logic nanodevices fea- 
ture high switching speeds, low power consumption, ex- 
cellent on-off current ratios and high scalability and in- 
tegrability. Electron transport in ultrasmall nanodevices 
approaching channel lengths of order ~ 10 nm [l], how- 
ever, turns out to be nonballistic and intricate [21 [3] 
due to intrinsic phenomena (such as thermal/shot noise 
caused by the discreteness of charge) and extrinsic per- 
turbations (such as scattering off unscreened trapped 
charges). These effects give rise to significant cur- 
rent fluctuations at high clock speeds and low voltages 
01 0, which are detrimental to efficient device opera- 
tion. Highly doped drain and source regions can fur- 
ther impact channel electrons, e.g., due to the build-up 
of mirror charges, suggesting that Coulomb interaction is 
a paramount ingredient in describing transport, dissipa- 
tion and equilibration in nanostructures [2j [31 [6] . A com- 
prehensive understanding of electron dynamics on small 
time and length scales therefore is desirable to improve 
the device performance. 

Various approaches to electron transport in nanode- 
vices have been taken so far, ranging from classical Monte 
Carlo and molecular dynamics methods to quantum 
non-equilibrium Green function calculations IH [5l [^ • 
Here we address the charge transport from an interme- 
diate, semiclassical perspective [9l [10], by solving the 
Schroedinger equation numerically for a pair of inter- 
acting electron wave packets propagating in a planar 
nanochannel. This approach interpolates between the 
classical, particle-like and the quantum, wave-like nature 
of electrons (Fig. fl]) and approximately preserves the 
discreteness of the transported electron charge over ap- 
propriate (not too short) time scales. Wave packet ap- 
proaches have been considered previously, but relied on 
strictly one-dimensional structures and/or on noninter- 
acting electrons, or were based on approximate schemes 
such as time-dependent Hartree-Fock theory [ll.-.16j . 

Below we study in detail the effect of Coulomb repul- 
sion and an external electric field on wave packet propa- 
gation in a two-dimensional nanochannel. A main finding 



FIG. 1. (Color online) Schematic illustration of classi- 
cal (particle-based), quantum (wave-based) and intermediate 
electronic transport regimes applicable to a nanochannel. For 
intermediate system sizes a semiclassical approach based on 
a wave packet picture may be used, which is further strength- 
ened by an external field and/or Coulomb interaction. (A 
indicates the electron mean free path.) 



is that Coulomb repulsion redistributes the charge den- 
sity to the channel walls, which makes electron transport 
more sensitive to perturbations at the interface. The 
presence of a uniform electric field, mimicking a chan- 
nel potential, furthermore has a stabilizing effect on the 
wave packets, reducing the spreading along the channel 
direction and hence backing the more localized, classical 
particle-like picture often used in full scale device simu- 
lations. 

Wave packet time evolution — In order to study the 
propagation of two interacting electrons in a nanochan- 
nel, we solve the time-dependent Schroedinger equa- 
tion ihdr\ip{r)) = II{t)\iI){t)) numerically using a split- 
operator (Suzuki- Trotter) approach [TTHIQ] . The total 
Hamiltonian II{t) = i^kin('?") + ^int includes a kinetic 
and a Coulomb interaction term, i.e.. 
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FIG. 2. (Color online) Comparison of the time evolution of two noninteracting and two interacting electron wave packets in 
a two-dimensional nanochannel in the presence of an electric field applied in x-direction. (a), (b) Snapshots of the electron 
density p(r;r) for [7 = eV and 2.0 eV, respectively, (c), (d) The corresponding partially integrated one-dimensional electron 
densities px{x; r) and py{y; r) for [/ = eV and U = 2.0 eV, respectively, (e) Total charge transported across the half-way line 
at xo — 23.5 nm as a function of time r. The initial sizes of the Gaussian wave packets are Gr — ^a — 4.5 nm, and the electric 
field amounts to ^ = 0.4 MV/cm; only the incoming, initially left-hand side wave packet possesses initially a finite momentum 
k= (0.37r/a,0). 
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respectively, where cj (cr) creates (annihilates) a spinless 
electron at site r. Here the kinetic Hamiltonian i^kin(T) 
describes nearest-neighbor hopping processes with hop- 
ping amplitude t and is time-dependent through a Peierls 
phase factor to account for a uniform electric field E = 
{E^ 0) in x-direction. Coupling the electrons to an elec- 
tric field in this way, the channel potential along the x- 
direction can be modelled without spurious discontinu- 
ities associated with finite size boundaries in a periodic 
channel along x (at the expense of the conservation of 
total energy, which becomes periodic in time), whereas 
open boundary conditions are imposed along y. The lat- 
tice spacing and the two lattice vectors of the underlying 
square lattice are denoted by a, a^^ and a^y, respectively. 
The term H-^^^ represents the long-range Coulomb inter- 
action with interaction strength U . 

We are interested in the evolution of the two parti- 
cle state iV^(r)) = (1/V^) Er,r' ^(^^ rV)4cMo). where 
^(r,r';r) stands for the antisymmetric two-particle 
real space wave function at time r, normalized to 
Xlr r' l^(^5^^'^)P = 1- ^s initial condition we model 
the probability distribution of the electrons by two (mov- 



ing) Gaussian wave packets centered at R and R' with 
momenta k and k' and width a^, 



*(r,r';0) = ^[$k(r - R),$k'(r' - R')], 



(3) 



where ^[. . .] ensures proper antisymmetrization and nor- 
malization and <l>k(i*) ex exp[— |rp/(2cr^) + ik • r]. Below 
we consider only scenarios where the right-hand side wave 
packet is initially at rest (k' = 0) while the left-hand side 
wave packet initially possesses a finite incoming crystal 
momentum centered at k = (0.37r/a, 0)] (see the snapshot 
sequences in Figs. |2]and[3|. 

For a realistic set of parameters applicable to nan- 
odevices such as Si based structures, we choose t = 1 
eV, which corresponds to an effective band mass of 
m* = 0.763 me in units of the free electron mass mg, 
and a = 0.5 nm for the lattice spacing. Since the channel 
lengths of nanodevices will reach only a few tens of nm in 
the near future, we further base the wave packet simula- 
tion on a lattice of size 47.5 x 15.5 nm^ (corresponding to 
96 X 32 lattice sites), representing a rather long channel 
to reduce finite size effects. 

Interacting electrons in a planar nanochannel — The 
typical time evolution of a wave packet pair in the pres- 
ence of a moderate finite electric field E = {)A MV/cm 
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FIG. 3. (Color online) Effect of an electric field applied in x-direction for a fixed Coulomb repulsion magnitude, (a), (b) 
Snapshots of the electron density />(r; r) for ^ = MV/cm and 0.4 MV/cm, respectively, (c), (d) The corresponding partially 
integrated one-dimensional electron densities Px{x;t) and pyiy'^r) for ^ = MV/cm and E — 0.4 MV/cm, respectively, (e) 
Total charge transported across the half-way line at xq — 23.5 nm as a function of time r. The initial wave packet sizes are 
given by fJr = 3a = 1.5 nm, while U — 2.0 eV; only the initially left-hand side wave packet is initially given a finite momentum 
k= (0.37r/a,0). 



and in the absence (U = 0) or presence {U = 2.0 
eV) of Coulomb interaction is summarized in Fig. |2l 
where a^ = 4.5 nm. The snapshots in Fig. |2] (a) 
and (c) for [7 = 0, which display the electron density 
/)(r; r) = 2^^, |^(r,r^;r)p and the corresponding one- 
dimensional electron densities px{x]r) = ^yp{Y]T) and 

PyiV'^^) = T^xPi^'^^) (^ = i^^y)) along the x- and y- 
direction, respectively, clearly imply that the evolution 
of the wave packets is strongly determined by the lat- 
eral confinement in y-direction and to a lesser extent 
by the Pauli exclusion principle (since we are consider- 
ing spinless electrons only). Due to the confinement in 
y-direction, the spreading of the wave packets and the 
subsequent reflection at the channel walls lead to inter- 
ference along y, e.g., causing a pronounced temporary 
maximum on the channel centre line at r = 16.46 fs. 
The partial reflection of the incoming wave packet off 
the initially static wave packet due to Pauli exclusion 
also causes an, albeit much weaker, interference pattern 
along the x-direction. Most of the incoming wave packet, 
however, passes by the initially static wave packet such 
that both wave packets remain essentially distinct and in- 
tact, propagating with different and increasing velocities. 
Under the present electric field, the peak-to-peak ampli- 
tude of the corresponding Bloch oscillation amounts to 



2so = At/{eE) = 100 nm, which substantially exceeds 
the channel length [20 . 

The effect of finite Coulomb repulsion {U = 2 eV) is 
clearly visible in the snapshot sequences of Fig. |2] (b) 
and (d). Due to the Coulomb interaction the incoming 
wave packet pushes the initially resting wave packet for- 
ward and enhances the transverse expansion along the 
y-direction. The latter results in a more pronounced in- 
terference pattern along y, with the overall probability 
weight of the wave packets pushed further toward the 
channel edges. The Coulomb interaction also slightly en- 
hances the charge transport along the longitudinal chan- 
nel direction, as can be seen from the total charge trans- 
ported across the line at xq = 23.5 nm [see Fig. [5] (e)]. 
The total transported charge Q{t) can be obtained by 
integrating the probability current density in x-direction 
over time and the channel width. The x-component of 
the current density operator is given by 



fij;T) 



ita , 



[e-'t^-cU^Cr-h.c.], (4) 



which can be derived from the appropriate continuity 
equation [21^, yielding 

Q{r) = (-e)^ r(^(T')|r(xo,y;T')|V(r'))dr'. (5) 



We furthermore note that the time-dependence of the 
total transported charge seems largely independent of 
the initial size, shape (e.g. quasi-one-dimensional Gaus- 
sians versus circular Gaussians) or positions (centered or 
slightly off-centered with respect to each other) of the 
wave packets. 

Wave packet stabilization in an electric field — The 
effect of a moderate electric field in the presence of fi- 
nite Coulomb repulsion ([/ = 2.0 eV) is displayed in Fig. 
[3J which compares the wave packet evolution foTE = 
and 0.4 MV/cm. Here we have chosen a pair of initially 
smaller wave packets with (Jr = 1.5 nm, which spread 
more quickly with time and generate a more heavily mod- 
ulated density distribution in the transverse channel di- 
rection [cf. panels (c) and (d) of Figs. |2]and[3]. The wave 
function nevertheless maintains a larger weight near the 
channel edges than in the centre, similar to the findings 
above. 

More strikingly, however, the overall expansion of the 
electron density in x-direction is drastically reduced in 
the presence of a finite electric field [cf. the snapshots 
in Fig. [3] (a-d) at time r = 26.33 fs], where a density 
maximum in x-direction clearly develops within the sim- 
ulation period. A similar electric field effect can also be 
observed in the case of larger initial wave packets. The 
electric field, expectedly, also leads to significant increase 
in the transported charge [Fig. [31(e)]. In the absence of 
an external field, the transported charge starts to satu- 
rate near unit charge |e| at about r = 20 fs, indicating 
that, despite the (weakly) current enhancing influence of 
the Coulomb interaction, largely only one electron wave 
packet reaches the right half of the channel during the 
simulation time. 

Discussion and conclusion — By modeling the initial 
electron probability distribution via two interacting wave 
packets, the present study takes a semiclassical approach 
to charge transport in a planar nanochannel [9l [10] . Our 
simulations show that the uniform spreading of the wave 
packets is limited along the transverse channel direction 
due to both interference effects and Coulomb repulsion. 
As a result, the electron density is strongly modulated 
along y-direction and propagates along the channel x- 
direction with a relatively large probability weight near 
(but not too close to) the channel edges. This density 
pattern can also be interpreted as the pair of electrons oc- 
cupying only the low lying transverse eigenmodes whose 
probability distribution is modified and pushed closer to 
the channel edges due to the Coulomb interaction. The 
details and the extent of such edge dominated transport 
can have significant implications for the performance of 
nanodevices, where interface roughness and charge trap- 
ping at interfaces are known to lead to substantial degra- 
dation of the electron mobility ^ ^5^ . 

Furthermore, we find that for typical nanochannel 
length and time scales an electric field has a stabiliz- 
ing effect on electron wave packets, inhibiting spreading 



along the channel direction. This behavior reinforces the 
more particle-like picture of electrons in a nanochannel, 
lending support to particle based classical or semiclassi- 
cal studies such as classical Monte-Carlo and/or molec- 
ular dynamics methods to solve the classical Boltzmann 
transport equation in the presence of interactions and 
impurities. The stabilizing effect of an electric field has 
also been observed for strictly one-dimensional interact- 
ing wave packets based on a time-dependent Hartree- 
Fock approach, which additionally revealed a bunching 
phenomenon of wave packets within the one-dimensional 
channel pT] , 

In conclusion, we have considered the effect of Coulomb 
repulsion and external field on the electron dynamics 
in a nanochannel, by exactly solving the Schroedinger 
equation for a pair of interacting and propagating elec- 
trons. This is in contrast to previous wave packet studies 
which mostly focussed on single electrons and/or one- 
dimensional nanostructures p!2HT6] . The present wave 
packet simulation may be extended in future studies to 
include, e.g. the electronic spin degree of freedom, small 
three-dimensional nanostructures, more than two inter- 
acting electrons and, though more challenging, the cou- 
pling to high density source and drain contacts. 

The authors thank Y. Tokura for useful discussions. 
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